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Abstract 

We investigate the nature of transition rates for the hadronization 
process of qq — ► htxtx as opposed to the transition rates for qq — ► -kit, 
within the Nambu-Jona-Lasinio model that has manifest chiral sym- 
metry. Feynman diagrams appropriate to this process can be classified 
according to the expansion in the inverse number of colors 1/N C . Two 
of these types of graphs are seen to be either "s-like" or "t-like" in 
nature. A further graph that contains both s-like and Hike elements, 
and which is denoted as si-like, is also present. To describe such a 
process with two incoming and three outgoing particles, it is necessary 
to extend the number of Mandelstam invariants. It is convenient to in- 
troduce seven such variables, of which only five are independent. The 
cross section for two incoming particles to three outgoing particles is 
then reexpressed in integral form in terms of these invariants. As a 
function of y/s, the final expression is reduced to an integral over the 
four remaining invariants. The limits of integration, which are now 
non-trivial, are also discussed. The transition rate for the explicit case 
of uu — ► 7r + 7T~7r°, is evaluated numerically, assuming non-chiral pions, 
m n = 135MeV. The rate for three pion production is found to be of 
the same order of magnitude as for two pion production, making this 
a non-negligible contribution to hadronization. 
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I. INTRODUCTION 



Nonequilibrium formulations of field theories are in a stage of development at 
present. Their study is essential for ultimately constructing a solid theoretical basis 
for the transport codes that are used according to various models for the simulation 
of heavy-ion collisions. One such attempt that is still in its early stages investigates 
the nonequilibrium formulation of a chiral theory, the Nambu-Jona-Lasinio (NJL) 
model From practical experience, one assumes that a collision term of the 

associated semi-classical limit will be Boltzmann-like in form, and this has been 
demonstrated explicitly pi. In particular, a full calculation that includes mesonic 
degrees of freedom must also contain not only quark- quark and quark- ant iquark 
scattering cross-sections on the partonic level, but also the hadronization processes 
qq —>■ MM' and their dissolution, where M represents a mesonic state. Cross- 
sections for these processes in SU(3) have already been calculated and are given 
elsewhere J3J. 

Processes such as those described in the last paragraph are binary in nature. 
On the other hand, there are numerous other processes that may occur and which 
increase the total number of particles. While binary collisions lead to thermal equi- 
librium, other multiparticle collisions may also be required to attain chemical equi- 
librium. In particular, the simplest possible further kind of collision would be one 
that increases the number of particles, such as that of qq —>■ MM'M" . Such a pro- 
cess should also rightly enter into the collision dynamics. If one had some knowledge 
of the magnitude of such processes, one would at least know whether multiparticle 
production in transport models could be expected to play an important role or not. 
In particular, one imagines that multiparticle hadronization processes such as the 
type just mentioned could be particularly important for the Goldstone bosons, the 
pions, since they are massless in the chiral limit, and are light in any event. Since 
pions are the dominant degree of freedom in the late phase of heavy-ion collisions, 
being produced copiously, it is essential to study their role. The purpose of this 
paper is thus to examine the cross-section for processes of the kind qq — > tttttt, 
where two particles enter the interaction region, while three are emitted. A direct 
numerical calculation is performed specifically for the case uu — > 7r + 7r~7r°. This spe- 
cific example has been chosen because the hadronization channel into two mesons 
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uu — > 7r + 7r has already been calculated and has been shown to dominate over the 
other possible existing two pseudoscalar emission channels that include uu — > 7r°7r°, 
K+K- and 7r°r/ [|]. 

The chiral model underpinning our calculation is the NJL model, in which the 
SU(2) version is used. While this model has several shortcomings, specifically that 
it does not confine, it gives an excellent description of the mesonic and baryonic 
sectors @. In using this model, one is able to obtain a self-consistent description 
within a chiral theory, which ultimately may be used in developing a self-consistent 
transport theory that in turn is based purely on this chiral Lagrangian [BL|5L|6|]. 

As a guide to determining the most relevant diagrams, the expansion in the 
inverse number of colors 1/N C is used |7|||. Within the same order of 1/N C , there 
are three types of diagrams, an s-like, a t-like, and a mixed so-called st-like graph. 
We find expressions for the invariant amplitudes in all three cases, and apply this to 
a calculation of the transition rates specifically for the process uu — > 7r + 7r~7r°. Since 
there are many terms that contribute to each channel, we investigate all channels 
first separately, and find that the st-like channel is at least an order of magnitude 
larger than the s-like and t-like graphs. That is to say, within the 1/N C ordering, 
there appears to be a subordering that dictates that one (or more) graphs are in 
fact dominant. In retrospect, one can give heuristic arguments why the particular 
st-like graph should in fact dominate. This particular graph contains the most 
divergent functions, since it contains a fermion loop with three propagators. The 
s-like graph by contrast contains a fermion loop with four propagators. In fact, such 
an ordering of terms within one specific class of diagrams that superficially are all 
of the same order in 1/N C rather according to the degree of divergence of each graph 
was proposed long ago in Ref. 0. The calculation in this paper thus serves also in 
part to validate a posteriori that this ordering correctly assesses the importance of 
each term. In Ref. |9[], convergent terms were regarded as being secondary, while 
the leading behavior was determined solely by the divergent terms. In a further 
calculation of the pion radius [[l(J , it has also been seen that the leading logarithmic 
contribution that is in accord with the well-known results of sigma models and chiral 
perturbation theory, is also due to the most divergent graph of the many available 
to the same order in the 1/N C expansion when calculated in the NJL model. This is 
so because this graph naturally corresponds to the pion loop diagram that has been 
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known to give rise to the chiral logarithm in so many other models or theories, see 



for example Refs. pTHT^I . Note that in Ref. ||10|| , however, no explicit calculations 
of the remaining terms was undertaken. 

In our numerical calculation of the transition rate that was carried out using 
a Monte Carlo procedure, we find that the st-like channel has the same order of 
magnitude for the production rate of uu — > 7r + 7r~7r° as is found for the dominant two 
particle production rate uu — > 7r + 7r~ Q. As such, it appears to give a nonnegligible 
contribution to the production cross-section for three pions and cannot be ignored 
in hadronization programs at low energies. 

An essential feature of the calculation is the construction of the cross-section 
for the process of two incoming to three outgoing states, a + b — > c + d + e. To 
this end, we introduce extended Mandelstam-like variables that will be denoted as 
s, t, u, v, w, x and y, and we reexpress the cross-section in terms of these variables. 
The invariant amplitude, which is given by the relevant Feynman graphs is most 
easily expressible in terms of these extended variables Q. Although this section of 
our work is performed for the case of incoming quarks and outgoing pions, and 
has thus the restriction pf = p\ = m 2 and v\ = v\ = v\ = m 2 , the kinematics 
are otherwise unrestricted so that the generalization to other arbitrary three body 
processes is obvious. 

This paper is organised as follows. In Section II, we evaluate the invariant 
amplitude for two arbitrary incoming quark flavors and outgoing mesonic states. In 
Section III, we discuss the kinematics of the process a + b — > c + d + e and write 
the cross-section for this process in terms of extended Mandelstam-like variables. 
In Section IV, we calculate the cross-section for uu — > 7r + 7r~7r°, and present our 
numerical results. We summarize and conclude in Section V. 

II. HADRONIZATION INTO THREE PIONS. 

In this section, we discuss the processes that lead to three pion production, 
governed by the SU(2) Nambu-Jona-Lasinio Lagrangian density 



1 Other choices are also possible [14]. The variables chosen here however, arise naturally 



when one examines the possible i-channel processes. 
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C = ${i?- m )^ + G[(^) 2 + (^7s^) 2 ] • (2.1) 

Here = uiq = mo is the current quark mass for u and d quarks and G is a 
dimensionful coupling strength, [G] = [MeV] -2 . There are several detailed reviews 
of this model in the literature, see for example Refs. 0. The reader is referred to 
these references for a derivation of the gap equation that determines the constituent 
quark mass 

d 4 p 



m = m + 2iGN c Nf J — ^trS(p) (2.2) 



(2*) 

with Slp)^ 1 —]f> — m as well as for obtaining the pseudoscalar and scalar meson 
masses via the condition 

l-2GYi s M k2 )\ m l, m l = 0, (2.3) 

where the scalar and pseudoscalar polarization functions are defined in the usual 

way as 

Tl PS (k ) = ^N c tr f (T lTf ) J (^i [p2 _ m2][(p + fc)2 _ m2] (2-4) 
tt / 7 2\ • a/- /• rf V (y* + m)(yH-ff + m) 

1 W ) - ll\ c J (2?r)4 [f _ [(p - fc)2 _ (2.b j 

with Ti an isospin operator T 3 = r 3 and T± = (n ± zt 2 )/a/2. 

Using the concept now that one can construct the Feynman diagrams for any 
process involving quarks, antiquarks and mesons, one can classify the graphs that 
contribute to hadronization. In the case that one has hadronization into two parti- 
cles, this is given in the standard fashion according to s channel, t channel and u 
channel graphs, as is shown in Fig. 1. Here s, t and u are the usual Mandelstam 
variables, s = (pi + £> 2 ) 2 , t = {p\ — p 3 ) 2 and u = (pi — pi) 2 . From energy-momentum 
conservation, p 1 + p 2 = Pz + Pi, one has the additional constraint 

s + t + u = 2m 2 + 2m\. (2.6) 

These amplitudes and their resulting cross-sections have been discussed in detail in 
Ref. ||. Note that the set of graphs is chosen to be leading in the 1/N C expan- 
sion Jr|J§. That the s, t and u channels are of the same order lies in the fact that 
meson propagators (in the pionic channel) are constructed as intermediate states 
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into quark-antiquark scattering amplitudes. These quark-antiquark scattering am- 
plitudes contribute as /m^ ~ l/iV c , as has been explicitly shown in Ref. ||, for 
example. A fermionic loop on the other hand brings a contribution of order ~ iV c . 

Generic diagrams that we can construct for two incoming quarks and three out- 
going mesons are shown in Fig. 2(a)-(c). We note that in direct analogy to the 
graphs of Fig. 1, we can order the hadronization graphs of Figs. 2(a) and 2(b) as 
being s-like or t-like. Figure 2(c) fits into neither category and we thus denote it as 
st-like. As has already been stated, the graphs of Fig. 2 are generic: the full set of 
diagrams is obtained for Figs. 2(a) and 2(b) by permuting the outgoing momenta. 
There are 3! = 6 ways of doing this. Whether all 6 diagrams in each case contribute 
to a process or not is then determined by the flavor factors of the individual in- 
coming quarks in question. In Fig. 2(c), four cases arise on noting that the vertex 
for two particle production can be attached either to the incoming (top or bottom) 
quark or antiquark line, and the sense of the line in the quark loop can be clockwise 
or anticlockwise. The single emanating meson line can also take on any of the final 
state momenta, leading to a further permutation of these three possibilities. Exact 
combinations again depend on the flavors of the vertices in question. 

For two incoming states to three outgoing states, it is useful to define the ex- 
tended Mandelstam-like variables, 



s = 


(pi 


+ P2? 








t = 


(pi 


-Ps) 2 








u = 


(pi 


-P,f 








V = 


(pi 


-P.Y 








w = 


(pi 


-P3- 


P,f = 


(Pa 


-P.? 


X = 


(pi 


~Pi- 


P5) 2 = 


(Pa 


-Ps) 2 


y = 


(pi 


~P3 ~ 


P,f = 


(Pa 


-P4) 2 , 



where the last relations follow from energy-momentum conservation, p% + p 2 = P3 + 
Pi + Pb- Using, in addition, the restriction to quarks for the incoming states, 

p\ = p\ = m 2 (2.8) 

and pions for the outgoing states, 
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2 2 2 2 

P 3 =P4=P5= 



(2.9) 



we find 



s + t + u + v = 3m 2 + 3m 2 



(2.10) 



and 



s + w + x + y = 3m 2 + 3m 2 



(2.11) 



in analogy to Eq. (|2.6|) . 

In what follows, in our analysis of the s-like and t-like graphs, the products P1P2, 
PiPs and so forth, occur. These are simply related to the extended Mandelstam 
variables via, e. g., 



We do not explicitly rewrite our expressions in terms of these variables: the nu- 
merical program makes a direct translation. In contrast to this, the t-like channel 
amplitudes are most simply written directly in terms of these variables and we thus 
use them explicitly from the start. 



In this section, we discuss the s-like channel diagram involving an intermediate 
scalar or pseudoscalar meson. One of the six possible ways of distributing the 
momenta is shown in Fig. 3. We shall assume that the final state always contains 
three pseudoscalar particles. In this case, if the intermediate meson is a scalar, the 
scattering amplitude, vanishes, 



This must be so on general grounds, since the parity of the final state is negative, in 
contrast to the positive parity of the intermediate state. This can also be validated 
by direct calculation. 




Pi -P3 = —=(t ~m 2 - ml). 



(2.12) 



A. s-like channel 



iM a (p 3 ,P4,P5) = 0. 



(2.13) 
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We thus require Fig. 3 calculated with the pion in the intermediate state. We 
examine this single diagram first with the momenta as indicated in the figure. Using 
an obvious notation, this is 

- iM v (p3,P4,Ps) = ^75?>i4 lC 2 1 _ 2GUp S (k 2 ) ^ 9 ^ A ^^- K ^ p ^' Ph ^ 

(2.14) 

where the four pion vertex is given as 
A 7T7T7T7r (p 3 ,p4,P5) = 

-Tr J -^^i^ 5 T k iS(q - k)i^ 5 T 5 iS(q - p 3 ~ p^i^T^iSiq - p 3 )i'j 5 T 3 iS(q). 

(2.15) 

In this expression, Tr = tr 7 tr/tr c sums over all degrees of freedom, and T, is an 
appropriate flavor operator. Color appears trivially as shown by the 5-function in 
Eq. fl2.14D . Inserting the single particle propagator leads to the form 



^77777777(^3,^4,^5) = -iV c ivj 345A:) x 

d^q tr 7 [(- 4+ jj + m) (£- p / 3 -p' 4 + m)(- 4+ j> 3 +m)(4+m)} 
(27r) 4 (q 2 — m 2 ) {(q — p 3 ) 2 — m 2 ] [(q — p 3 — p 4 ) 2 — m 2 ] [(q — k) 2 — m 2 ] 

(2.16) 

and the flavor factor is abbreviated as jvj 345 ^ = tr ^(TkT^T^T^). The flavor factor to 
be associated with the complete s-like channel is then 

/f 45 ) = jvfivf 45fc) , (2.17) 

where N9 = qTjq arises from the choice of incoming quark and antiquark. Directly 
performing the spinor trace leads to the expression 

^77777777(^3,^4,^5) = -2N c Nf Abk) [(-k 2 p 3 p 4 - kp A ml)I 4 (k,p 3 +p A ,p 3 ,0) 
+2(kp 3 )pJ 4 (k,p 3 +p 4 ,P3,0) +p 3 pj 3 {p 3 +P4,P3,0) 
+ (kp 4 - k 2 )I 3 (p 3 + p 4 , k, 0) + (kp 4 - p 3 p 4 - ml)I 3 (p 3 ,p 3 + p A , k) 
+h(P3+P4,0) + I 2 (k,p 3 ), (2.18) 

where 
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d A q 1 

InK ...,a n )= I (27r)4 [( 9 _ ai )2_ m 2]...[( 9 _ an )2_ m 2] ( 2 ' 19 ) 



and J4 is denned as 

h( ai , . . . , = J (27r )4 [ (g _ fll )2 _ m 2] . . . [(g _ Q4 )2 _ m 2] " ( 2 - 2 °) 

This latter integral, in combination with an external momentum can be reduced to 
a linear combination of the I n . These functions are listed in Appendix A. 

The diagram of Fig. 3 represents one of the six possibilities of ordering the 
external momenta. Another of these is obtained by reversing the arrows on the 
internal loop of the vertex, as shown in Fig. 4. Translation of this Feynman diagram 
gives 

^7r7T7T7r(P5,P4,P3) = 

-Tr J -^^i75T k iS(-q)i^ 5 T 3 iS(^^ 

(2.21) 

Some algebra is required to evaluate the traces. The spinor trace can be shown to 



be equivalent to that required for Eq. (|2.16|) . In addition, the flavor traces are equal 



Af 43fc ) = Af 45fc ), (2.22) 

so that 

Ar7r7r7r(P5,P4,P3) = ^tttttttt G»3 , Pi , P5 ) • (2.23) 

This is a general feature when reversing the loop momenta for the diagrams required 
here. Hence the combination of six terms for the s-like channel reduces to a cal- 
culation of 2 times three terms. The remaining two "crossed" s-like graphs for the 
vertices are shown in Fig. 5. The associated vertex functions are A n7Tnv (p4,p 5 , p 3 ) 
and A vvn7T [p^, 

P4,P3), which can be obtained from Eq.( |2.18| ) on suitably interchang- 
ing the arguments. 

One may finally construct the scattering amplitude for the s-like channel as 
-iM s ~^ c = -2i[M w (p 3 , Pi ,p 5 ) + M w ( Pi ,p 5 ,p 3 ) + M n (p 5 ,p 3 , Pi )}, (2.24) 
where 



9 



2G 

- iM w (p a ,p b ,p c ) = v 2 nbUi 8 C1 >C3 1 _ 2 ^, k2 . gl m N^ ) A mrmr (p a ,p b ,p c ). (2.25) 

Since IIp5(/c 2 ) and can also be expressed through the integrals I n , 

- tU PS (k 2 ) = 47V C (27 1 (0) - A; 2 / 2 (0, -A;)) (2.26) 
g^ q ^-AiN c I 2 (0,-k) : (2.27) 

a direct evaluation of this quantity depends only on the evaluation of the integrals 
I n , see Appendix A. 

For the purpose of obtaining a scattering cross-section, the absolute value of the 
scattering amplitude is averaged over initial states and summed over the final ones. 
Since Z) S1 S2 |^75Wi| 2 = 2s, it follows that 



2G 2 g$ qq (Nf) 2 s\A 7T7T7T7T (p a ,p b ,p c ) 
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"v 'jf^m (2 - 28) 



from which |.Mf-~ like | 2 can be constructed. 



B. t-like channel 



The six possible t-like channel graphs are shown in Fig. 6, together with our la- 
belling of each process. As an example, we calculate one of these diagrams explicitly. 
Considering the t — w channel depicted in Fig. 6a, one finds 

- iM {t ' w) = v 2 ig^ qq i^T^iS{pi - p 3 - Pi)ig^ qq i^T A iS(pi - p^ig^i^T^Sc^, 

(2.29) 

and therefore the absolute value of the scattering amplitude averaged over initial 
states and summed over the final ones follows as 

1 _f2 6 AT 

4N Z £j 2 (t — m z ) z {w — m z ) z 

•(fa- i>z + ^)75MiMi75(^i- ^3 + "1)75(^1- ^3- A + 771)75]. 

(2.30) 

The trace over flavor matrices has been abbreviated as ft jW = ^{T^T^T^). After 
some tedious algebra for evaluating the spinor trace, one arrives at the relatively 
simple form 
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f 2 o 6 



qq 



2N r 



-1 



t — m? 

2 



x — u 

1+ n I " 



rrit 



w — m 2 ) {t — m 2 ) 2 {w — m 2 y 



m 



y — m 



u — m 



(2.31) 



(t — m 2 )(w — m?) \ t — m 2 w — m? 

in which the invariants t, u, x, y and w occur. Completely analogous calculations 
lead to the expressions for the remaining five channels, lA^'^l 2 , |.M(' U ' 2/ )| 2 , \Ai^' y ">\ 2 , 
\JH( U > W )\ 2 and lA^^'^^l 2 . These can be obtained directly from Eq.( |2.31[ ) by making an 
appropriate variable substitution. To obtain [.M*"'^! 2 , one makes the replacement 



t — > u — > v 
w — > x — > y — > w 



(2.32) 



in Eq.(|2.3l|), that follows on examining the momentum dependence of Fig. 6b in 



comparison with that of Fig. 6a. Similarly, the expression for lA^'^l 2 follows from 
Eq. fl2.31|) on substituting 



t — > v u — > t 
w — > y — > x — ► w, 



while |A^(*'^| 2 is obtained on writing 



(2.33) 



t — > t x — > x 
u — > v w <-> y, 



(2.34) 



\M^\ 2 from Eq.( [2~3lD on writing 



t <-> M 



and lA^^I 2 from Eq. ([23l| ) on writing 



(2.35) 



t — > f w — > w 
w <-» x y — > y. 



(2.36) 



All expressions have relatively simple analytic forms. One notices that in the chiral 
limit, m,r = 0, they simplify dramatically. 
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C. si-like channel 



The st-like channel of Fig. 2(c) is evaluated by separating the vertex for two 
meson production from that of single meson production. This diagram is thus sub- 
divided into two pieces, the first containing the quark-qntiquark scattering amplitude 
that is mediated by an exchanged meson connected with a quark loop as indicated 
in Fig. 7. Implicit in this diagram is that both the clockwise and counterclockwise 
senses in the fermion loop are to be evaluated. Since we will regard pions the final 
state, it follows from parity conservation that the intermediate mesonic state must 
be a scalar. In this case, this is then taken to be the a meson. Since the two flavor 
NJL model is considered here, the two pions that emerge from this vertex must be 
charge neutral. Translating Fig. 7, one has 

A awn (p) = iD s ip 2 )(ig wqq ) 2 {-2N c )tT f T 3 T 4 
d A q 



( 2n y tT ^5iS(q)i-f 5 iS(q + p 3 )iS(q - p 4 )]. 

(2.37) 

The factor 2 arises from considering both senses of the fermion loop, and again tr 7 is 
the trace on Dirac indices. In this expression, Ds(p 2 ) is the quark-quark scattering 
amplitude in the scalar a channel, i.e. 

D ^ = 1-2W - < 2 ' 38) 

In principle, the exchange of the a meson within the model as required here, involves 
a complex function, since the a is a resonance. It is, however, only weakly unbound, 
with m 2 ~ 4m 2 + m 2 , and consequently here only principle values of the integrals 
I n entering into the evaluation of the irreducible polarization Us(p 2 ) is used. 

By making the substitution q — > — q in the second integral of Eq. (|2.37|) and 



performing the trace, one arrives at the following expression for this function, 

A awn (p) = 16miV^ 2 (?9J D s (p 2 )[/ 2 (p3,-p 4 )+p3.p 4 J3(0,p3,-P4)]. (2.39) 

This graph is then embedded into the hadronization graph of Fig. 2(c). Attaching 
the vertex A a7T7T to both the incoming quark line, as was drawn originally in Fig. 2(c) 
plus the incoming antiquark line, which is also a required configuration, leads to two 
terms, 
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{P5—P2)-m 2 

= A^ st) + M 2 st) . (2.40) 

The averaged transition matrix element then follows as 

1 1 



tu — m 2 f — m 2 



x {m 2 (4m - s) + (v — m — m^)(w — m — m 2 )}. (2.41) 

Since p^ is the momentum of the pion line that is attached singly, we use the notation 
M-^iPh) to denote this fact. 



D. Full scattering amplitude 



The full scattering amplitude is comprised of all s-like, t-like and st-like contri- 
butions, i.e. 

M = 2M n (p 3 ,p 4 ,p 5 ) + 2M 7T (p 4 ,p 5 ,p 3 ) + 2M 7T (p 5 ,p 3 ,p 4 ) 

+ + M {u ' x) + M (v ' y) + M {t,y) + M iu,w) + M {v ' x) 

+ M {st) (p 3 ) + M {st \ Pi ) + M (st) (p 5 ). (2.42) 

As such, the cross-section requires not only the averaged terms that were given 
in Eqs.( [2.28D and fl2.31| ) - ( |2.3(j| ), but also mixed terms that arise on building the 
modulus. There are many of these and we have listed some of them in Appendix B. 



III. THE SCATTERING CROSS-SECTION FOR qq tttttt IN TERMS OF 
EXTENDED MANDELSTAM-LIKE VARIABLES 



A. Cross section 



The notion of relativistic invariance is an important calculational guide in eval- 
uating matrix elements. As such, it is particularly useful to express the integrated 
cross section in terms of relativistic invariants only. To do this in the case when two 
particles are incident and two are exiting is a standard textbook exercise, see for 
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example ||15|| . However, when three particles are involved in the exit channel, as is 
the case here, this becomes somewhat more complex, and we are required to derive 
a formal expression. The only restriction that occurs in our derivation is that the 
incident particles with momenta p\ and p 2 refer to particles with the same mass, 
i.e. p\ = p\ = m 2 , while the exiting particles with momenta p%, p 4 , and p^ also are 
bound by the constraint p\ = p\ = p\ = m 2 . The generalization to arbitrary mass 
particles in describing the reaction of a + b—>-c + d + e is obvious. 
Starting from the definition of the cross section 



a 



(27r) 4 £( 4 )Q 3 +p 4 +p b -pi -p 2 ) \M(p l7 p 2 ,P3,P4,P5)\ 2 
\v re i\2Ei2E2 
d 3 p 3 d 3 p 4 d 3 p 5 
X (2vr) 3 2E 3 (27r) 3 2E 4 (27r) 3 2E 5 1 j 

and the definition of the extended Mandelstam-like variables of Eq. (|2.7|) , we use the 
identity 

Jdp 5 e(p 5 )5(pl-ml)2pl = l (3.2) 
and perform the p$ integration, giving 



f \M(s,t,u,v,w,x,y)\ 2 d 3 p 3 d 3 p 4 2 2 

a = J (2^m^E-2E2 ^m 5iiPl+P2 - p3 ' Pi) 

X e( P $+P2-P3-P4)- (3-3) 

From energy conservation, the argument of the G function can be evaluated to be 
Pi + P2 ~ Ps ~ P4 = Ei + -^2 — -^3 — E± = E 5 > 0, so that the function takes the 
value 1. The argument of the 5 function can be reexpressed in terms of the extended 
Mandelstam-like variables either as 

(Pi + Pi ~ P3 ~ Pi) 2 — m \ = s + w + x + y — 3m 2 — 3m 2 (3.4) 

or as 

(Pi + P2 - P3 - Pa) 2 -ml = s + t + u + v- 3m 2 - 3m 2 , (3.5) 



and the fact that the 8 function contributes only when the right hand side of Eqs. ( |3.4| ) 
or (|3.5| ) is zero leads to the energy momentum relations already stated in Eqs.( [2.i0| ) 



or 
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As a consequence of their definitions in Eq. (|2.7|) , the variables s, t, u and w are 
unaffected by the p§ integration. On the other hand, v , x and y are altered, 

v = ( Pl - P5 ) 2 _> ( P2 -p 3 - p 4 ) 2 
x = (pi-p4- p 5 ) 2 -> (P2 - P3) 2 

2/ = (pi -P3 -P5) 2 -> (pa -P4) 2 - (3.6) 

The next task lies in transforming the element of integration and determining the 
appropriate Jacobian. To this end, we denote 6^4 as being the angle between the 
vectors p\ and p 3 or p 4 . Since we work in the cm. reference frame, p\ and p 2 lie 
in opposite directions so that the angle between p 2 and p 3 or p 4 is (ir — 6> 34 ). In 
addition, \pi\ = \p 2 \, so that the extended Mandelstam-like variables of Eq.fl2.7p take 
on the forms 

t = m 2 + ml - 2E 1 E 3 + 2\p\\ \p 3 \ cos9 3 (3.7a) 
u = m 2 + ml - 2E X E 4 + 2\px\\p 4 \ cos8 4 (3.7b) 
v = m 2 + 2m 2 n ~ 2E i E 3 ~ 2E i E 4 + 2E 3 E 4 
-2|pi||p 3 | cos 6» 3 - 2|pi||p 4 | cos6> 4 

— 2\p 3 \ [p 4 |(sin 6 3 sin 6 4 cos{ip 4 — ip 3 ) + cos 9 3 cos 84) (3.7c) 
w = m 2 + 2m 2 n ~ 2E i E 3 ~ 2E i E 4 + 2E 3 E 4 
+2\p i .\\p 3 \ cos 6» 3 + 2|pi||p 4 | cos6> 4 

— 2|p 3 | |p 4 |(sin 6 3 sin 6 4 cos(</9 4 — ip 3 ) + cos 9 3 cos 84) (3.7d) 
x = m 2 + ml - 2E X E 3 - 2\p x \\p 3 \ cos9 3 (3.7e) 
y = m 2 + ml - 2E X E 4 - 2|pi||p 4 | cos# 4 . (3.7f) 

For conciseness, we have used the abbreviation E 34 = |p 3i4 | 2 + ml, but draw the 
attention of the reader to the fact that we consider |p 3 4 | together with the angles 
as variables, and not £3,4, althougth this would be equivalent. By contrast, E\ and 
E 2 are fixed by the centre of mass condition, E\ — E% — y/s/2. Conversely the 
inverse transformation can be obtained. By adding and subtracting Eq. (|3.7a|) from 
Eq.( [3.7e|) , one can easily find the determining equations for |p 3 | and 6 3 . A similar 
manipulation of Eqs.( |3.7b] ) and ( |3.7tp yields those for \p 4 \ and 9 4 . Having done this, 
Eqs. (|3.7c|) or (|3.7d|) can be easily solved for cos(y9 4 — ^3). The exact form of this 
transformation is 
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t — X 

t - x = 4|pi||p 3 | cos 6» 3 |p 3 1 cos 6> 3 = — — , (3.8a) 

Mpi\ 

u - y = 4|pi||p 4 | cos6» 4 =^> I&|cos0 4 = -^-p (3.8b) 

9 9 2m 2 + 2m 2 — t — x 

t + x = 2m 2 + 2ml ~ ^E 3 E 3 = , (3.8c) 

AE± 

o 9 2m 2 + 2m 2 — u — y 

u + y = 2m 2 + 2ml ~ 4 ^i^4 => E A = — £ (3.8d) 

while the expression for cos <y2 34 is 

t + u- w-m 2 + 2E 3 E 4 - 2|p 3 ||p 4 | cos0 3 cos0 4 , _ 

cosy? 34 = — - — ■ — . (3.9) 



Returning to the integrand in Eq.( |3.3| ), one sees that it depends only on the 
extended Mandelstam-like variables, that in turn can be written in terms of | ^0*3 1 , 
\Pa\, 3 , 64 and the difference (p 34 = <p± — ip 3 . Thus 

J dip^dipi = 2n J d(pu- (3.10) 

There are five independent integration variables, \p 3 \, \p 4 \, 8 3 , #4 and (^34. On the 
other hand, there are seven extended Mandelstam variables. As can be seen from 
Eqs.( |3.7c"| ) and ( |3.7cl| ), both v and w are functions of the variable <£> 34 ; we arbitrarily 
choose to discard v. In keeping with convention, the other unrelated variable is 
s, the square of the center of mass energy of the incoming particles. We thus are 
required to relate the multidimensional volume element in the variables |p 3 |, |p 4 |, 
# 3 , 6 '4 and <y9 34 to the extended set t, u, w, x and y, which follows as 

2\J\ d\p 3 \ d\pi\ dOzdO&dtfM = dt du dw dx dy, (3-H) 

where the factor 2 has been introduced together with the reduction of the integration 
region of 9934 from [0, 27r] to [0, 7r]. After some calculation, one finds the Jacobian of 
the transformation to be 

|J| = 128 |pi| 2 |p 3 | 3 |p 4 | 3 sin 2 # 3 sin 2 ^ 4 ^r|sin^34|. (3.12) 



Using the information from Eqs.( |3.4| ) and ( 3.12 ), the cross section takes the form 



. \M'(s,t,u,w,x,y)\ 2 f 2 2N 

a= l (2^4^ S{s + w + x + y-ton -3m,, 



\P3 


2 

sm 


#3 




2 sin 


9 4 E 3 E 4 dt du dw dx dy 


2E 3 2E 4 256\p 1 \ 2 \p 3 \ 3 \p 4 


3 sin 2 # 3 sin 2 8 4 E\ 


sin 9934 1 



X om oTn <ircl - 191- 131- I3„:_2/1 „:„2 /l m2 I' (3.13) 

(3.14) 
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which can be written in a simpler form as 
1 



a{s) = 64(2^ J dt dU dW dX dy N(s(s - 4m 2 ))3/=» *' % W ' X,V) ^ 

x 6(s + w + x + y — 3m 2 — 3m 2 ), (3.15) 



where use has been made of the relationship 4|p*i| 2 = s — 4m 2 . In both Eqs.( p.l^ ) 



and Q3.15 )) the notation M' has been introduced to remind us that the variable v in 
the calculated matrix elements is to be replaced by (— s — t — u + 3m 2 + 3m 2 ) and 
the factor 



iV = 2[p 3 |[p4| sin ^3 sin ^yl — cos 2 (^34 (3.16) 

has been abbreviated in the denominator. Clearly iV must also be expressed in 
terms of the extended Mandelstam-like invariants. One finds the rather unwieldy 
expression 



N = N(s, t, u, w, x, y) 
-1 



[(2m 2 + 2m 2 ) (t - x - u + y) + lux - 2ty\ 



As(s — Am 2 ) 

2 f (t - x) 2 + (u - yf (2m 2 + 2m 2 - t - x) 2 + (2m 2 + 2m 2 - u - yf 



\ s - 4m 2 
— {t + u — w — m 2 ) 2 + 4m^. 

+ {t + u — w — m 2 ) 



(I ■<■)[„ y)^ 

s — Am? 



(2m 2 + 2m 2 — t — x){2m 2 + 2m 2 — u — y) 
s 

(3.17) 



In the chiral limit, this function simplifies somewhat. 

As a final step, the integration on w can be performed and we obtain our final 
expression for the cross section: 



f x 1 f A+ A r , \M"(s,t } u } x,y)\ 2 

a{s) = J dt dU dX dy {*{* - Am 2 ))V 2 N(s, t, u, x, y) ' (3 ' 18) 

which involves a fourfold integration. The notation M" serves here to remind one 
that w has been integrated out, and is to be replaced by (3m 2 + 3m 2 — s — x — y) 
in these matrix elements. Using an obvious notation, the function 

N(s, t, u, x, y) = N(s, t,u,w = 3m 2 + 3m 2 — s — x — y,x,y) (3.19) 
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can be obtained from Eq.( 3.17|) . This result generalizes the well-known form for the 



cross section for two incoming particles to two outgoing particles, 



a{ [JM^M dt (320) 

v ; J 167rs(s-4m 2 ) v ; 

in which Ai'(s, t) denotes matrix elements M'(s, t) = A4(s, t, u = 2m 2 + 2m 2 T — s — t), 
for qq — > irir say. 



B. Limits 

In reactions of the type a + b — > c + d, for which the cross section of the form 
Eq.( |3.20D holds, the determination of the limits of integration of the single variable 
t is simple: it is always related to one angular variable by virtue of its definition, 
and the minimum and maximum values of this angular variable lead to a minimum 
and maximum value of t respectively. For the case that is implicit in Eq. (|3.20|) , a 
simple analytic form for t max and t m i n can be found. One has 



t max = m 2 + ml-- + Vs- 4m 2 ^/- - ml (3.21) 



and 



tmin = m 2 + ml - - - y/s- Am 2 J --ml (3.22) 

for the case of different particles and 

tmin = m 2 + ml- I (3.23) 

when there are identical particles in the final state. 

By contrast, in our case where there are three outgoing particles in the final 
state, we have four integration variables t, u, x and y, whose ranges have to be 
determined, and which are laid fixed by the maximum and minimum values that 
these variables can take according to Eqs. ( |3.7a| )-( |3~TIf ) which are in turn set by 
the appropriate minimal and maximal values of the angular and energy variables 
that occur therein. For this reason, it is difficult to give simple expressions for these 
limits and we thus start by giving their ranges. Obvious constraints on the particle 
energies and angles involved are 
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E 3 > m n 
E 5 > m n 



cos 2 e 3 < 1 
cos 2 e 4 < 1 

COS 2 </?34 < 1 



(3.24a) 
(3.24b) 
(3.24c) 
(3.24d) 
(3.24e) 
(3.24f) 



These constraints on the angles must be translated into conditions for t, u, x, and y, 
while an expression for cos<£> 34 can be obtained on inverting Eq.( |3.7d|) for w. Now, 
introducing the abbreviations 



t = t — m 2 — m 2 



u = u — m 2 — mt 



x = x — m 2 — m 2 



y = y — m 2 — m 2 



(3.25a) 
(3.25b) 
(3.25c) 
(3.25d) 



we can reexpress Eqs.( p.24"a| ) and (|3.24b|) as 



t — x — 2m n \/s > 



(3.26) 



— u — y — 2m^s > 0, (3.27) 

which follows on using the equations ( |3.8cj ) and ( |3.8d| ). The equation ( |3.24cj ), on 
the other hand, becomes 



s - E 3 - E A - > 0, (3.28) 

on using the conservation of energy condition E3 + E4 + E5 = in the cm. system. 
E 3 and E4 can be eliminated from Eq.( |3.28| ) using ( |3.8c| ) and ( |3.8d| ). Then it reads 



r- t + x u + y , . 

^ + vf + v?- m - £0 - (3 ' 29) 

This expression can be further utilized to set an upper bound on either t + x or 
u + y. Let us regard that for t + x. Eliminating u + y using Eq. ( |3.27| ) , Eq. ( |3.29| ) 
leads to the bound 



19 



t + x > -2s + Ay/s: 



(3.30) 



We next examine the angular condition Eq.( |3.245 ). In terms of t and x, or t and 
x, this reads 

2 

I <1, (3.31) 



t — x 



\ 4 |Pl||P3 

where Eq. (|3.8a|) was used. By eliminating 4 1 j?\ | 2 = s — Am 2 and \f>z\ 2 = E\ 



ml 



with the aid of Eq.( |3.8c( ), this condition becomes 

(t- xfs <{s- Am 2 )((t + x) 2 - Asm*) 



(3.32) 



In an analogous fashion, the angular condition Eq. ( |3.24e|) can be rearranged to give 
a condition on u — y: 



u — y) s < (s — Am )((u + y) — Asm^ 



(3.33) 



The final angular constraint, Eq.( 3.24f ) is rewritten in terms of the Mandelstam-like 
variables on considering Eq.( |3.9| ), replacing w by 3m 2 + 3m 2 — s — x — y and the 
appropriate energies from Eqs.( |3.8c|) and (|3.8d| ). This then takes the form 

, 2 



(s + t + u + x + y + m\ + aj/s — flS/ii) 



A{El - m\){E\ - ml) - 26(£f - m 2 )/5 - 2/3 {E^ - m 2 )/5 + [3 2 5 2 /s 2 
where the abbreviations 



=?<1, (3-34) 



a 



7 



t + x 
u + y 



(3 
5- 



t — x 
u - y 



(3.35) 
(3.36) 



y/2' V2 

have been introduced and s = s — Am 2 . 

The integration region for the variables t, u x and y must be determined from 
the relations flpp , (ggg ), (|3^0|) , O and From (gg§) , it follows 

immediately that 



a < 



25772, 



(3.37) 



while Eq. ( 13.32 ), considering the equal to sign, describes a hyperbola in the variables 
a and /3 with defining equation 
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The hyperbola has semi-axes a = y2sm% and b = J2sm%, and asymptotes (3 = 
±ba/a. In view of Eg. ([37371 ), only the negative hyperbola must be considered, which 
is depicted in Fig. 8. In order to determine whether the integration region lies within 
or without this function, it suffices to choose a point in the plane. Choosing (3 = 
and a < —a leads to the relations 

* = -*7T = -;75 (3 ' 39) 



so that 



2(1 + ^E 3 ) = V2p = (3.40) 



(s - Am 2 )(El - ml) - (t + ^E 3 ) 2 > 0, (3.41) 

i.e. the condition Eq. (|3.32j) is fulfilled. The integration region thus lies within the 
hyperbola. The bound onHi that was found in Eq. (|3.30| ) can be translated to the 
variable a and gives rise to the vertical line 

a > -V2s + A^ m7T (3.42) 

shown in Fig. 7 with dots and which cuts off the integration regime to the left. This 
defines a clearly delineated area F\ in the a — j3 plane for the integration region. 
Note that one could transform the conditions on a and j3 back to conditions on t 
and x. This results in a hyperbola for these functions that is obtained on rotation of 
the a and /3 axes by 7r/4. It is not necessary to find this form explicitly, because the 
checks that are done to evaluate the integral will determine the range of integration 
numerically using the conditions on a and (3. 

In an analogous fashion, it is possible to construct an area F<i in the variables 
7 — 5 that are related to u and y that is delineated by the analogous condition 

7 > -V2s + 4^m n (3.43) 

plus a hyperbola in the 7 — 5 plane. Taken together, these two areas F\ x F 2 span a 
superset of points in t-u-x-y space that constitute the region of integration, subject 
to the condition Eq. ( |3.34j ) being fulfilled. 
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IV. NUMERICAL CALCULATIONS 



In this section, we examine the reaction uu — ► tt + 7t~7t°. This particular process 
has been selected, because the analogous two particle final state hadronization pro- 
cess uu — > 7r + 7r~ is the dominant channel Q for uu —>■ MM', where M, M' are 
mesonic members of the SU/(3) multiplet. In order to make a sensible comparison, 
we calculate the transition rate for this process, which is given in terms of the cross 
section a(s) as 

u(s) = \v re i\a(s), (4.1) 
where v re i is the relative velocity of the incoming quarks, and which is given as 



KU=2^1-^ (4.2) 

in the centre of mass frame. Our numerical calculations are based on a Monte Carlo 
procedure for evaluating the fourfold integral in the expression Eq. (|3.18 ) for the 



cross section. In implementing this, it proves useful to alter the integration region. 
As already mentioned, this consists of surfaces F% and F2 in the a-(3 and j-S plane 
respectively and is bounded by the further condition Eq. fl3.34p . First we extend 
the surfaces F\ and F2 to the enlarged regions W\ and W2 that contain them, but 
which have simpler boundaries. The triangular region bounded by the asymptotes 
of the hyperbolae together with the vertical restriction on the value of the ordinate 
as indicated in Fig. 8 for the (3-a plane is chosen as the extension of F% to W±, and 
a similar construction is made for W2 in the j-8 plane. This has the advantage that 
random points can be chosen to be equally distributed in this triangular area for 
the Monte Carlo procedure. On actually performing the numerical integration, the 
caveat is imposed that the function to be integrated takes the required functional 
value in Fx, F2 but which is zero otherwise, together with the restriction Eq.( |3.34j ). 
The delineation of the integration regions for a and (3 is given via 

C^max ^ <Z 
-Pmax < (3 < Pmax (4.3) 

with a max = 2/a/2 - y/s/ v^m^ and (3 max = a max V§/ y/s. 
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In order to evaluate the transition rate, the individual contributions to A4 that 
occur in Eq. (|2.42|) must be specified. These are in turn dictated by the traces of the 
flavor factors, which sometimes vanish. For the specific combination of 7r + 7r~7r° in 
the final state, the s-like channel contribution contains all six possible graphs, 

M s - likc = 2M 7T (P3,P4,Ps) + 2M 7T (Pi,Ps,P3) + 2M 7r (p 5 ,p 3 ,p,), (4.4) 
while in the t-like channel, only three possible graphs contribute, 

M t-hke = M (t,w) + M (v, y ) + M {t,y)_ (4.5) 

In the st-like channel, there are only four possibilities. Here the single outgoing 
meson must be a 7T°, while the loop contains the 7r + and 7r~ outgoing in the final 
state. The two loop directions plus the attachment to either the incoming quark or 
antiquark line comprise these possibilites. They are fully summarized in the function 

M st - likc = M (st) ( P5 ). (4.6) 

Rather than evaluate the complete transition rate that is constructed from A4 = 
j^s-hkc _|_ ^t-hke _|_ st -hko ^ we eV aluate first the transition rate in each particular 

channel, in order to illustrate the magnitude that each type of channel contributes. 
We comment briefly on our numerical procedure. In the s-like channel, the pseu- 
doscalar polarization function Ups(k 2 ) is required from Eq.( p,14|) , while for the 
st-like channel, the scalar polarization function Ils(A; 2 ) enters as in Eq.( |2.39| ), cor- 



responding to the exchange of it or a mesons respectively. At zero temperature, 
the case studied here, one has that at k 2 = m 2 , the tt is strongly bound, while at 
k 2 = m 2 , the a is weakly unbound, since m 2 ~ 4m 2 + m 2 . Thus the a already has a 
small width, representing the unphysical decay into two quarks. At higher values of 
k 2 , this unphysical width grows, and at large values of k 2 , the pion can also develop 
a width. We have chosen to ignore such unphysical decays into quark constituents. 
This is done primarily for reasons of numerical simplicity - only the principle values 
of the integrals I n that enter into Ups(k 2 ) and Us{k 2 ) are considered - and this 
allows one to focus on the role of the triangular vertex graphs, which are the phys- 
ical cause of the decays into the mesons. This procedure simply assumes that the 
intermediate mesonic states are stable with regard to decay into their quark con- 
stituents, which is not unreasonable. Note that the imaginary parts of the I n that 
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enter into the irreducible polarization functions commence always at the threshold 
k 2 = 4m 2 . A detailed discussion of these functions giving their analytic structure, 
together with graphic illustrations as a function of k% and k 2 can be found in Ref. 



For our numerical results, the parameters A = 851 MeV, GA 2 = 2.87, and 
m = 5.2 MeV have been used. These lead to the values f n = 93 MeV, m n = 135 
MeV and (ipip) = (—250 MeV) 3 . Calculated values of the quark mass yield m = 265 
MeV and the pion quark coupling as being g nqq = 2.85. 

Figures 9,10, and 11 display the results for the s-like, t-like, and st-like channels 
individually. One sees that the st-like channel provides the largest yield, being 
substantially larger than the s-like channel (an order of magnitude) and the t-like 
channel (two orders of magnitude). Thus, one sees that within the set of graphs that 
are of the same order in the 1/N C expansion in fact have a very different relevance. 
It is important to try to understand this, and we can do so on a heuristic level: 
one notes that both the s and st like channels contain a quark loop. The fewer the 
number of quark lines occurring in this loop, the more divergent, or less convergent 
this diagram will be. Graphs such as that occurring for the s-like channel contain 
a quark loop with four internal quark lines in contrast to the graphs of the st-like 
channel, which has a quark loop with three internal quark lines. One may thus 
expect that the st-like channel will dominate over the s channel. Our findings and 
this argument are in line with that of Ref. JJ], in which all convergent terms are 
dropped and only the leading divergent or least convergent ones are retained. A 



similar argument was given in Ref. [ID| for ignoring several 0(1/N C ) terms in the 
higher order corrections to the pion radius, in which the same argument was made 
on heuristic grounds. This numerical calculation now validates such heuristic claims. 
Thus this ordering plays an additional role in reordering the set of terms that occur 
beyond the 1/N C expansion. 

Obviously, by examining the graphs of Figs. 9 to 11, the full cross-section should 
be well approximated by that of the si-channel alone. This is given in Fig. 12, and 
includes the cross terms. We also compare our result with the transition rates for 
the processes uu — ► 7r + 7r~ and uu —>■ 7r°7r° that proceed via the diagrams of Fig. 1. 
Here the process uu — > 7r + 7r~ is the dominant one. We find, for the range of validity 
of the model, say ^/s < lGeV, that the hadronization rate of u and u into three 
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pions, which is also driven by a two pion hadronization graph, is of the same order 
of magnitude as uu into two pions, despite the fact that the phase space is altered. 
This result indicates that transport or other models that include hadronization as 
a two meson level, require hadronization into at least three pions also. 



V. SUMMARY AND CONCLUSIONS 



In this paper, we have studied the three meson hadronization process qq — > 
MM' M" within the two flavor Nambu-Jona-Lasinio model. Our motivation stems 
from the fact that pions are light, almost massless particles and can thus be produced 
copiously. We have calculated transition rates for this process in order to assess the 
importance of this process in comparison with the standard binary processes that are 
usually input into transport models. Since the NJL model is currently under study 
as a transport model for a chiral theory PJ5|JT6||, such questions can be answered 
within the context of the model itself. It has thus been our aim to investigate the 
transition rate for this process, in comparison with similar transition rates for two 
meson production q — > MM' for a specific initial state q and q being u and u. 

In our analytic discussion, we have found that there are three types of diagrams 
that contribute in the lowest level 1/N C expansion that is appropriate to this model. 
Our numerical calculation however indicates that an additional criterion for deciding 
on which graphs dominate the set of diagrams in a particular 1/N C class is the 
number of internal fermion lines when a fermion loop diagram is present. This 
criterion is in accord with the heuristic methods of Ref. 0, and which was used 
successfully in identifying the logarithmic contribution to the pion radius [IU|, and 
thereby recovering the known results of sigma models, other chiral models and chiral 
perturbation theory (CHPT) for this quantity. We have thus explicitly been able in 
this paper, to demonstrate that this crucial heuristic argument given in [10] is in fact 
valid, in particular in the context of a different calculation. In view of the fact that 
the dominance of the most divergent graph gives the expected result for the pion 
radius, our result here may not simply be a model artifact. We have at present no 
way of knowing how this result would compare with that calculated within another 



framework, say a confinement model [13] or within QCD itself, where the mechanism 
giving rise to the leading term may be different. We avoid any speculation here. 
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We conclude by stating that the numerical calculation shows that our results 
lie within the same order of magnitude as the hadronization rates for two pions, 
thereby indicating that such processes must necessarily be taken into account in 
hadronization at low energies. 
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APPENDIX A: PAULI-VILLARS REGULARIZATION FOR THE 

INTEGRALS I N 



The evaluation of I n {a\ . . . , a n ) as defined in Eq. (|2.19|) is done using the Pauli- 
Villars regularization scheme. The reader is referred to one of the many papers 



involving such integrals for I\ and I 2 , see for example ||8| JI7| , [18| . We simply quote 
the final results, 



(4tt) 2 

and I 2 (jp, 0) = J 2 (p) is 

Up) = 



' [ m 2 \n(l-( n A2 ,J ^+2A 2 lnfl+ A ~' 



m? + A 2 ' 



77T 



+ A 2 



(Al) 



2 ,1 

]T Cj / dz\n{m 2 - z(l - z)p 2 ). 



(A2) 



W 2 U 

where Co = 1, C\ — 1 and C2 = —2 is a standard set of coefficients fulfilling 
J2j Cj = J2j Cjm'j = 0. Here m 2 = m 2 + a 2 A 2 , with ao = 0, a\ = 2 and a 2 = 1. 
Reference |19| gives 



h(p,k) 



^2Cj / dyy \ dx[ 



S7r "0 J ° 



a' 



2(a 2 



\2J' 



(A3) 



with a = y(l — x)p + (1 — y)k and bj = y(l — x)p 2 + (1 — y)k 2 — m 2 , and 

3i 2 f 1 f 1 f 1 1 
Ii(p,k,l) = — — A , - Cj / dzz dy dxx 2 [— -rrr, tttt + 



(An) 



j=0 



JO 



"J- 
a 2 — bj 



4(a 2 -^) 2 6(a 2 -^) 3 



, (A4) 
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with a = zx(l — y)p + x(l — z)k + (1 — x)l and bj = zx(l — y)p 2 +x(l — z)k 2 + (1 — 
x)l 2 — m|. An analytical expression for I 3 can be found in fll9| , while ^4 has been 
evaluated numerically. 

APPENDIX B: MIXED TERMS IN THE SCATTERING AMPLITUDE 

There are many mixed terms that occur between the various channels. Here 
we list only those that directly enter into the calculation of the transition rate for 

MM — > 7T + 7T~7T°. 



1. Mixed terms between t-like graphs 



In this subsection, we list the cross terms that arise between t-like graphs in con- 
structing the averaged matrix element, and which are required by virtue of Eq.( 
There are three such terms. The first of these is 



J\4 (t,w) M(v,y)* 



ft,wfv,y9-jr 



4N r 



1 1 (u — m 2 ) + (x — m 2 ) 

y — m 2 t — m 2 (v — m 2 )(w — m 2 ) 



+ 



(t — m 2 )(x — m 2 ) 



+ 



(u — m 2 )(y — m 2 ) 



(v — m 2 )(w — m 2 )(y — m?) (t — m 2 )(v — m 2 )(w — m 2 ) 



+ m 
+ 



+ 



+ 



(v — m 2 )(w — m 2 ) (v — m 2 )(y — m?) (t — m 2 )(v — m 2 ) 
-2(t — m 2 ) — (x — m 2 ) y — x — 2{u — m 2 



{v — m 2 ){w — m 2 )(y — m 2 ) (t — m 2 ){v — m 2 ){w — m 2 ) 



+ ml 



{t — m 2 ){v — m 2 ){w — m 2 )(y — m 2 ) 
the second that is required is 



(Bl) 



ft,wft,y9TT 



+ 



l it — m 



+ 



4N C [ t — m 2 (t — m 2 )(w — m?) 
(v — m 2 ) (x — m 2 ) 



(t — m 2 )(y — m 2 ) (w — m 2 )(y — m 2 ) 
1 1 



{t — m 2 ){w — m 2 ) (t — m 2 ){y — m 2 ) (t — m 2 )' 



+ 



(x — m 



(w — m 2 )(y — m 2 ) (t — m 2 )(w — m 2 )(y — m 2 ) 
2 1 



— m 



w (t — m?) 2 {w — m 2 )(y — m?) 



(B2) 
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and the last one is 



fv,yft,y9w 



(w — m 2 ) 



4:N r 



+ 



[x — m 



y — m? (t — m 2 )(y — m?) 
(u — m 2 ) 



(v — m 2 )(y — m?) (t — m 2 )(v — m 2 ) 



+ 



+ 



2u — w — x 



(y — m 2 ) 2 (t — m 2 )(v — m?) (t — m 2 )(v — m 2 )(y — m?) 
2 1 



— m 



11 (t — m 2 ) (v — m 2 ) (y — m 2 ) 2 



(B3) 



2. Mixed terms between s-like graphs 

Mixed terms between s-like graphs do not take a simple form in terms of the 
extended Mandelstam-like variables. Generally, one has 

M n (p a: p b: p c )Ml(p' a: p' b ,p' c ) 

7T7r7r7r (p aiPbi Pc) A n7Tn7T 

- 8G 97Tqq (N f ) 5 C1C2 \i - 2GU PS (k 2 )\ 2 ' (B4) 

where {p a ,Pb,Pc) or (p' a ,p' b ,p' c ) represent the combinations (p3,P4,Ps), (PiiPsiPz) or 
(p 5 ,p 3 ,p 4 ). Averaging over incoming states and summing over final states leads to 



(B5) 



M w (p a , p b , p c )Ml(p' a , p' b , p' c ) 

2G 2 g 6 wgq (N}) 2 4sA nnnn ( Pa , Pb , p^A^Jp^ p' b , pQ 
N c \l-2GU PS (k 2 )\ 2 

which is evaluated directly numerically. 

3. Mixed terms between s-like and t-like graphs 

Here three such products are required. The first of these is the averaged ma- 
trix element of M. n (Pa,Pb,Pc) together with its complex conjugate, which can be 
constructed as 



M^{p a ,p b ,p c )M^ w > 



1 MGgl qq NjA 

7r7T7r7r (p aiPbiPc) 9tt qqft,w 

= N~csus2 V2ll5Ul 1 - 2GU PS {k 2 ) (t - m 2 ){w - m 2 ) 
xm 1 7 5 T 5 (^i- i> z + m)7 5 r 4 (^i- j^ 3 — ^ 4 + m)^T z v 2 - (B6) 
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Upon taking the traces and forming the complex conjugate, one finds the expression 



M^(p a ,p b ,p c )M^ w > + M^M*{p a ,Pb,Pc 

G9®qqft,wNj I A 7T7T7Tn (p a ,Pb,p c ) \ 2 2 ._, 



2Ke [ oATrr — j—^- ] \(y — m — m v )(w — m — m v ) 

(B7) 



N c (t - w?) (w-m 2 ) \ 1 - 2GU PS (k 2 
(t — m 2 — m n )[(2(w — m 2 ) + (x — m 2 — 2m 2 K ) — [u — m 2 — m w )(y — m 2 — 2m^)} 



and correspondingly 



M*(pa,Pb,Pc)Mt v *>+M<- v >tiM*(p a ,Pb,Pc) = 

G 9*mf">y N f 2 Re ( A ™^ a,Pfe,Pc ^ Uu - m 2 - m 2 )(v - m 2 - m 2 ) 

N c (t - m 2 )(w - m 2 ) K€ { 1 - 2GU PS (k 2 ) J i[U JU/ J 

+ (v — m — m 7r )[(2(y — m ) + (w — m — 2m n )\ — [t — m — m n )(x — m —2m n )}, 

(B8) 

and 



M v (pa,Pb,Pc)M<fi>V> +M^M* n {p a ,p b , Pc ) = 

G ^ggfv,y N f 2R ( A w ^(p a ,p b ,p c ) \ 2 _ 2 2 _ 2x 

+(i — m — m 7r )[(2(?/ — m ) + (x — m — 2m^)\ — (v — m — m n )(w — m —2m n )} 

(B9) 

which ends this subsection. 



4. Mixed terms between si-like and i-like graphs 

Using the generic labels Pa, Pb and pc to denote the outgoing pions in the t-like 
channel, one is required to form the cross term between the expression 

M = {ig n ) fv{p 2 )i^- - -275- -i^ 5 u[pi), BIO 

[Pc — P2) — m (Pi — Pa) — Tnr 

where / is a flavor factor, and the expression for M st ~ hkc [ n Eq.( [2.40| ). The averaged 
cross product using the first term of Eq. ( |2.40| ) gives 



1 Wnqq^cnnrf 



4N C (w — m 2 ) [(pi — pa) 2 — rn 2 } [(pc — P2) 2 — m 2 } 
tr 7 [(^i + m)7 5 (^i- $ A + m)7 5 (^ c - i> 2 + m)7 5 (^ 2 - m)'y 5 (p' 5 - i> 2 + m)]. (Bll) 
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The spinor trace can be performed, so that 



(at) -like »^t— like* _ y ^11 



Mr m 



m aa A <r™fm 1 



N c (w — m?) [(pi — pa) 2 — m?} [(pc — P2) 2 — m 2 ] 
x [-{Pa ■ Pc)Pb ■ (pi + Pa) + (pa ■ Ps)(pi - P2) • Pc - (pc • PsW ■ (pi - £2)]- 

(B12) 



Carrying out this procedure with the second term of Eq.( [2.40D leads to precisely 



the same spinor trace, so that the expression differs only from the above in that 
uj-m 2 is replaced by v-m 2 , in building A^^'^-M*- 1 * * and .M( s *)- likc .A/l*- like * = 

When the incoming states are u and u, and the outgoing states 7r + 7r~7r°, three 
t-like graphs can be constructed in which the outgoing momenta are permuted. 
These cases are (a) Pa = P5 = Ptt ? Pb = P3 = Pn+ and Pc = Pi = Pw- with / = 2; 
(b) p A = p 3 = P7T+, Pb = P5 = Ptt° and p c = p 4 = Ptt- with / = -2, and (c) 
Pa = P3 = Ptt+5 Pa = P4 = Ptt- and pc = P5 = Pn° with / = 2. Combining these 
three contributions and moving to the extended Mandelstam variables, one finds the 
form 



W(ai)-Uke A^t-like* — ^ m 9nqq A 



+ 



y; — m 2 f — m 2 



f — m 2 )((x — u — v — m 2 ) — m^x — y — v — m 2 ) 

X : 



k> — m 2 ]\v — m 2 ] 



(w — m 2 )(v — y) — (v — m 2 )(t — m 2 ) — m 2 (w — t — u + m 2 ) 



[t — m 2 ] [y — m 
(w — m 2 )(y — t — v + m 2 ) — m 2 (u — v — t 



m 



2\ 



+ ± ^ f . ^ % '-] (B13) 

in total. The mixed term between the (si)-like and s-like graphs does not have a 
simple form, and must be constructed from Eq.( [2.40P and Eq.( [2.28| ). 
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FIGURE CAPTIONS. 



FIG.l: s, t and u channel graphs available for the hadronization of qq — > MM', to 
lowest order in the 1/N C expansion. 

FIG. 2: s-like, t-like and st-like channels available for the hadronization of qq — > 
MM'M" , to lowest order in the 1/N C expansion. 

FIG. 3: Exchange graph in the s-like channel, for qq — > MM'M". 

FIG. 4: Vertex associated with the s-like channel hadronization graph. 

FIG. 5: The vertex 7r — > -r-r-k cross channels obtained by permutations of the external 
meson indices. 

FIG. 6: Six available t-like channels. 

FIG. 7: Two particle production from an intermediate state. 

FIG. 8: Region of integration in the a-j3 plane. The solid line indicates the negative 
hyperbola that is determined by Eq.( [3.38| ). The dotted vertical line is the boundary 
given by Eq.( |3.42j ). The asymptotes of the hyperbola are indicated by the dashed 
line. The surface Fx is bounded to the left of the solid line and the right of the 
dotted vertical line. W\ is the triangle bounded to the right by the dotted vertical 
line and to the left by the dashed lines. 

FIG. 9: Transition rate to plotted as a function of yfs for the s-like channel. 
FIG. 10: Transition rate u plotted as a function of \fs for the i-like channel. 
FIG.ll: Transition rate u plotted as a function of y/s for the st-like channel. 
FIG. 12: Transition rate u> plotted as a function of \/s. 
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